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This article presents an algorithm for estimating the signal- to- noise ratio (SNR) 
of signals that contain data on a downconverted suppressed carrier or fhe first 
harmonic of a square-wave subcarrier. This algorithm can be used to determine the 
performance of the full-spectrum combiner for the Galileo S-band (2.2- to 2.3-GHz) 
mission by measuring the input and output symbol SNR. A performance analysis 
of the algorithm shows that the estimator can estimate the complex symbol SNR 
using 10,000 symbols at a true symbol SNR of —5 dB with a mean of —4.9985 dB 
and a standard deviation of 0.2454 dB, and these analytical results are checked 
by simulations of 100 runs with a mean of -5.06 dB and a standard deviation of 
0.2506 dB. 


I. Introduction 

Current plans are for the performance of the full- 
spectrum recorder (FSR) and full-spectrum combiner 
(FSC) for the Galileo S-band (2.2- to 2.3-GHz) mission to 
be measured by running symbol-error-rate (SER) curves 
on the demodulated data. To do so, one needs to generate 
test signals at various points of the system, and the test 
signals must be consistent from one point to the other, 
which may be difficult. The measurement can only be 
done at the output of the demodulator, which means that 
if the demodulator is malfunctioning, other modules (e.g., 
the FSR or FSC) cannot be tested. Also, since test signals 
are needed, online testing is impossible. 

To overcome the disadvantages of measuring the per- 
formance from the SER curves, one would like to directly 
measure the symbol signal-to-noise ratio (SNR) at various 


points of the system. This means that the symbol SNR has 
to be estimated through complex symbols. By complex 
symbols, we mean the real symbols at an offset frequency 
with in-phase (I) and quadrature (Q) components. This 
article presents a complex symbol SNR estimator that es- 
timates the symbol SNR of data with I and Q components 
of a carrier and/or the first harmonic of a square-wave 
subcarrier. This SNR estimator can be used to measure 
the performance of the FSC as well as that of the complex 
symbol combiner for the Galileo S-band mission. 

The technique of the complex symbol SNR estimator is 
similar to that of the split-symbol moments estimator [1,2], 
except with complex symbols rather than real ones. The 
idea behind splitting a real symbol into two halves and 
correlating them (originally suggested by Larry Howard 
[2]), is to get the signal power by correlating two halves 
of a symbol, where the signal parts of the two halves of 
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the same symbol are correlated but the noises on the two 
halves of the symbol are uncorrelated. 

The same idea is applied here to the complex symbols, 
with the assumption that the offset frequency is known. 
The idea is to correlate the first half of a symbol with the 
complex conjugate of the second half of the same symbol, 
and then to take the real part to get the signal power multi- 
plied by a factor that is a function of the offset frequency. 
This factor can be divided out if the offset frequency is 
known. The same step is repeated for N symbols, and the 
results are averaged to get a better estimate of the signal 
power. 

The signal power plus the noise power can be obtained 
by averaging the magnitudes squared of the sums of the 
samples in a complex symbol over N symbols. The esti- 
mate of the SNR ratio is then readily obtained. In the 
next section, a brief description of the complex symbol 
SNR estimator is given. 

Some practical problems are not considered here. For 
instance, the frequency offset may not be known and may 
not be a constant. Also, if the symbol synchronization 
is off, say by a sample, the performance of the estimator 
will be affected, and in the future, the effect needs to be 
quantified. 

II. Brief Description of the Estimator 

For given samples of I and Q components of baseband 
data with an offset frequency lj 0 , we wish to estimate the 
symbol SNR. The I and Q components of the / th sample 
in the Mh symbol are expressed as follows: 

YYI 

Vhk ~ cos (“olT, + <}> 0 ) + n, lk (1) 

™ S 

yQ,„ = Jf d k sin ( l>JT s + <f> 0 ) + n Qtk (2) 

where 

/ = 0 , N s — 1 
k = 1 , N 

Here N s is the number of samples per symbol and is as- 
sumed to be an even integer, T s is the sample period, and 
dk = ±1 is the kth symbol. The noise samples in the I 
and Q channels are assumed to be independent with zero 
mean and equal variance <t^/N 8 . The true symbol SNR is 


snr = r ( 3 ) 

Presented here is an algorithm to estimate this symbol 
SNR, First, add all samples in the first half of a symbol 
for I and Q separately, 

N s / 2—1 

Y aIk = Y, W* (4) 

1=0 

N s / 2-1 

Y a Qk = Y ( 5 ) 

1=0 

Repeat the first step for the second half of the symbol 

Ns-\ 

Y Ph = Y V, ' k ( 6 ) 

l=N s /2 

N a - 1 

Y pQ k = Y y Q> k ( 7 ) 

l=N s /2 

Multiply the I component of the first half of a symbol by 
that of the second half; do the same for the Q components; 
and add the results. This is .equivalent to taking the real 
part of the product of the first half and the complex con- 
jugate of the second half of a symbol. Repeat the same 
procedure for N symbols, and average the results to obtain 
the parameter m p 

1 N 

m P = TV X] ( Y <*h Y Ph + Y aQk Y 0Qk) 

k=l 

= ( 8 ) 

Jt=l 

where 

Y ak = Y aIk +jY aQk (9) 

Y/) k = Yp Ik + jY/3 Qk (10) 


233 



Take the magnitude squared of the complex sum of the 
two symbol halves, repeat for N symbols, and average the 
results to obtain m ss 


= vX>«. +y <>. 


(ii) 


k= 1 


E {SNR*} ~ g (rn p ,m,,) + - 


( 2 ® 2 9 n / \ d 2 9 2 d 2 9 

* (»K +2c ° v + 


(14) 


Mathematical expressions for the complex symbol halves 
Y ak and Yfi> and the parameters m p and m ss are given in 
Appendix A in terms of the underlying model parameters 
of Eqs. (1) and (2). Finally, use m p and m ss to obtain the 
estimated symbol SNR 


where and cr^ ss are the variances of m p and m SS) 
respectively, and cov (m py m ss ) is their covariance. The 
variance of the estimate, SNR*, can be approximated by 
[3, p.212] 


SNR" = 

4 m p sine 2 (u 0 T s /2) / [sine 2 (w 0 T sy /4) cos (u 0 T 3y / 2)] 
m ss - m p (2 + 2/ cos (u> 0 T sy /2)) 

(12) 




<7 


2 

m>si j 


dg dg 
dm p dm. 


cov (m p ,m,,) 


(15) 


where w 0 is the offset angular frequency, and Tsy is the 
symbol duration, 


T —NT 

± sy — 


The estimate, SNR*, defined in Eq. (12) is shown in 
the next section to have an expected value equal to the 
true symbol SNR defined in Eq. (3), plus some bias that 
decreases with the number of symbols averaged. The vari- 
ance of SNR* is also derived in the next section as a func- 
tion of the number of symbols averaged and the true sym- 
bol SNR. 


In Eqs. (14) and (15), all the partial derivatives and the 
covariance are evaluated at (fn py 7n ss ). Expressions for the 
mean and variance of m p) the mean and variance of m 5J , 
the covariance of m p and m ssy and the partial derivatives 
are derived in Appendices B through E, respectively. Sub- 
stituting all the terms from the appendices into Eqs. (14) 
and (15), the mean and variance of SNR* are obtained, 


E {SNR"} = SNR-f 


1 

Ncos 2 (z/ 2) 


x 


^2 +SNR 



III. Mean and Variance of Symbol SNR 
Estimates 


r 2 

+ SNR 2 -J-(1 - cos (z)) 


(16) 


In Eq. (12), the estimate of the complex symbol SNR, 
SNR*, is a function of the random variables m p and m ss . 
Let g denote that function 


SNR* - g(m py m ss ) (13) 


and 


2 

^SNR* 


ATcos 2 (z/2) L C? + SNR C? +SNR2 


Assume that g(m py m ss ) is “smooth” in the vicinity of the x 

point (rn p ,m a5 ), where m p and m S3 are the means of m p 
and m ss , respectively. The expectation of the estimate, 

SNR*, can then be approximated by [3, p. 212] (17) 


(*-§«.«) snr 3 ^ ( i- 


cos ( 2 )) 
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where 


Z L0 O T S y 


and 


sine (u 0 T S y/A) 
sine (u> 0 T s /2) 


IV. Analysis Verification 

Two methods are used to verify the analytical results 
of the mean and the variance in the complex symbol SNR 
estimation: 

(1) Compare the numerical results of the analysis to sim- 
ulation results. 

(2) Set the offset frequency to zero and compare the 
performance of the complex symbol SNR estimator 
and that of the real symbol SNR estimator. 

A. Simulation Results 

Simulation results for the mean and variance of SNR* 
are shown in Figs. (1) and (2), and are compared to the 
theoretical results of Eqs. (16) and (17). The analytical 
and simulation results are also compared in Table 1. In all 
cases, the offset frequency, f 0 — u> 0 /(27r), is set at 20 Hz. 
The simulation results are based on 100 runs each, and the 
simulated means and variances are obtained by averaging 
over the 100 runs. 

B. Special Case 

A special case is when the offset frequency u> 0 = 0. In 
this case, the performance of the complex symbol SNR 
estimator should reduce to that of the real symbol SNR 
estimator given in [1] and [5], with N samples of pure noise 
and N samples of the signal plus noise. 

For this special case, the mean and variance of the com- 
plex symbol SNR estimates reduce to 

E {SNR*} = SNR + -^r(l + SNR) (18) 

oJ NR . = ^ (2 + 4SNR + SNR 2 ) (19) 

The analysis in [1], which assumed a constant signal level 
for all the samples, is easily extended to cover the case of a 


signal in only half the samples. The real symbol SNR and 
the real estimate, f?{SNR*}, in Eqs. (17a) and (17b) of [1] 
are simply replaced by values SNR and SNR* , respectively, 
representing averages over the 2 N real samples. These av- 
erage values are one-half the corresponding complex SNR 
and complex estimate SNR* defined here. With these cor- 
respondences, the mean of the real estimate in [1] reduces 
to 


£?{SNR*} = SNR + jy~y(l + SNR) (20) 


and the variance of the real estimate is 

j2 _( n y i 


X 


(SNR 2 + 2SNR + 1) + — (2SNR + 1) 


( 21 ) 


For a large A, Eq. (20) is equivalent to 

£{SNR*} — SNR + i(l + SNR) (22) 


and Eq. (21) approximates to 

4sr. = jf( SNR 2 + 4SNR + 2) (23) 

Comparing Eq. (18) with Eq. (22), and Eq. (19) with 
Eq. (23), it is clear that the special case of the complex 
symbol SNR estimator has the same performance as the 
real symbol SNR with N signal-plus-noise samples and N 
pure noise samples. 


V. Conclusions 

A complex symbol signal-to-noise ratio estimator is pre- 
sented in this article. This estimator modifies the split- 
symbol moments estimator for estimating the real symbol 
SNR [2] in order to accommodate complex symbols, with 
the assumption that the offset frequency is known. This 
estimator can be used to measure the performance of the 
full spectrum combiner as well as the complex symbol com- 
biner for the Galileo S-band mission. 
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Table 1. Simulation results of the complex symbol SNR estimation. 


Mean of SNR‘ , dB a of SNR* , dB 

SNR, dB N Runs 



Theory 

Simulation 

Theory 

Simulation 




-2 

-1.9968 

-2.0039 

0.2956 

0.2834 

2,500 

10 

100 

-2 

-1.9984 

-2.0136 

0.2112 

0.1937 

5,000 

10 

100 

-5 

-4.9985 

-5.0589 

0.2454 

0.2506 

10,000 

10 

100 

-8 

-7.9986 

-8.0437 

0.3056 

0.2787 

20,000 

10 

100 






Appendix A 

Sums of the Halves of a Symbol and Definitions of m p and m ss 


The sum of the first half of the fcth symbol is 


N,/ 2-1 

^~d k e- 

N s 

= E 

1=0 


- 171 Hkf 

J $o \ ~ 

- N, dkt 

1 




+ n c 


N s - 1 

nnp k = ^2 ( n /m + J n Qik) ( A ' 4 ) 

l = N s/2 


where j = a/ — 1 . 

Let s ak and sp k denote the signal parts of Y ak and Yp k , 
respectively. Then Y ak and Yp k can be expressed as 


= ^md k e^ 1/4WoT3 y 


l/aup ,.t, + ri sing jucTiy/A) 
sine (w 0 T 3 /2) 


and 


+ 


(A-5) 


+ n at 


(A- 1) 


Y P* = S P* + n 0y (A’ 6 ) 


and the sum of the second half of the Hh symbol is 


Ns - 1 


+ *•> + »* 


l=N, 


s/2 


= HL dk e^° T ° N3 ' 2 + 

JV S 


£o) 


1 _ e i{uoTsN,lt) 
1 _ e j« 0 T, 


+ n Pk 


_I m J e j( 3 /4^oT Jy - l/2w 0 Ts + 4„) SlnC O^y/ 4 ) 

2 * sine (u; 0 T a /2) 


The product of the sum of the first half of the Arth sym- 
bol and the conjugate of the sum of the second half of the 
same symbol is 

Yct k Ypk — S <*k S t3 k + S ^k n l3 k + n O( k s p k + n a k n p k 
= 1 m 2 c ,T-u>„T,„/2) Sinc 2 (u; o r, y /4) 

4 sinc 2 (Lj o r 5 /2) 

T s <y k T ri p k T n ot k s p k T na k rip k (A-7) 


+ n Pk 


(A-2) 


Taking the real part of the product, one obtains 


where the noise terms, n ak and np k , are complex Gaussian 
noises, with each of the components (real and imaginary) 
having a variance of cr^/2. In terms of the in-phase and 
quadrature noise samples in Eqs. (1) and (2), the complex 
noise terms, n ak and np k , are 





sine 2 (u 0 T, y /4) 
sinc 2 (w 0 T s /2) 

s } k + n ak n *t3 k ] 


(A-8) 


and 


Ns/ 2-1 

n ak = E ( n lik+j n Q>k) 

1=0 


(A-3) 


Define m p as 


m p 




Jb=i 


(A-9) 


23d 


The sum of the whole k th symbol is 


n-k ~ nj k + jnQ 


(A-12) 


n = E %-d k e^° T ’ 1 + M + n k 
1 = 0 $ 


where s/. fc , SQ k and n/ k , nq k denote the I and Q com- 
ponents of the signal and the noise in the k th symbol. 
Furthermore, Y* can be expressed as 


m 1 _ pju 0 TsNs 

= ™d k e>*- \ e +n k 


N, 


\ — (>jw 0 T s 


Yk = s Ik + js Qk + n Ik + jn Qk 


(A-13) 


md k A y ! 2 “° T ’y - + M si ” C ^°V 2) + „* 

sine {u 0 T s /2) 


— s k + 


(A-10) 


The sum squared of the kth symbol is 


,, r .2 o sinc“ (w 0 T ai/ /2) . ,2 , rtcvxr 

1**1 = m 2 . 2 V y / + K| +2#{n*«;j 

sine (a/ 0 T,/2) 


where s* is the signal part of the whole symbol and r\k is 
a complex Gaussian noise with a variance of in each 
component (real and imaginary), and 

Sk = s a k T 

Tik = n ak + ^/? fc 


where 


= m 2 C| + [n a- | 2 + 2${n*$J} 


,y 2 _ Smc2 ( UoT,y/2) 

to 


sinc"(u-’ 0 Tj/2) 


(A- 14) 


Note that the signal and noise parts of a whole complex 
symbol, s k and n kl can also be expressed in terms of their 
in-phase and quadrature components, that is, 


s k = sr k + jSQ k 


(A-11] 


Define m, 3 as 


= E i y * 


k=[ 


(A- 1 5) 
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Appendix B 

Evaluation of the Mean and Variance of m p 


All the expectations in the following are conditional on 
knowing the symbol synchronization. 


Because the {Uk} are independent, the variance of m p 
can be computed as [4, p. 352] 


k = 1 


1 N 
= N? Y 


a U k 


k=] 


(B-2) 


1. 2 sine 2 ( uj 0 T sy /4 ) 

4 sine 2 (uj 0 T s / 2) 


cos 


L 


Now let subscripts I and Q denote the in-phase and 
quadrature components of the signals. The second mo- 
ment of Uk can be expressed as 


1 N 

+ m22 + n “-4„ +n Qk n}J} 


k = l 


E i U k) - -7T E { s Ph + s ll k + s PQ k + s aQ k } 


1 2 sine 2 (lj„T s y /4) 

— 771 

4 sine 2 (u> 0 T a /2) 


cos 


UoTsy \ 


(B-l) 


(j 2 

+ -y + E{(sai k spi k + s aQk SfiQ k ) } 


where 

Uk = 


ro 2 sine 2 (cj 0 T, y/4) al <j\ 

2 sinc 2 (w 0 T s /2) 2 2 


= + n ak ){s*p k + n ieiJ} 


+ E {(s a i k spr k + s a Q k sp Qk ) 2 } 


and s ak and sp k denote the signal part of Y ak and Yp k , 
respectively, as in Eqs. (A-5) and (A-6). The last step 
in Eq. (B-l) follows from the independence of the noises 
with zero mean in the two halves of each symbol. Note 
that T sy = N S T S . 


Combining Eqs. (B-l), (B-3), and (B-2) yields 


JL f ^ -rn^ S ^ nC (^o^y/4) 2 a n \ 

N \4 Sine 2 (u> 0 T s / 2) n 2 ) 


(B-3) 


(B-4) 
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Appendix C 

Evaluation of the Mean and Variance of m„ 

ss 


The mean of m ss is 

E{m„} = E j^£l y *! 2 } 

= m 2 C 2 2 +^E^K | 2 + |E Emn k s' k }} 

V k=l V k=l 

= m 2 C 2 + 2<r 2 n (C-l) 

To find the variance of m SSi let 

V k =\Y k \ 2 

= K + 

= ( sj k + Ti[ k ) 2 + ( SQ k + nq k ) 2 (C-2) 

where Yj k , Yq h , s/ fc , sq* , n Ik , and ng fc are defined in Ap- 
pendix A. Because the {14} are independent, the variance 
of m ss can be expressed in terms of the variance of 14 [4, 
p. 352], 

E 0 Vi, (C- 3 ) 

fc=l 


4 fc = W 2 }-(£{^}) 2 ( c - 4 ) 

where 

E 2 {1/} = £ 2 {m JS } 

= (m 2 C| + 2<t“) 2 (C-5) 

Expanding I 4 2 , take the expectation of I 4 2 ; apply 
E{n) k j - E{riq k } = <r 2 , E{n\ } = £{n|J = 0, and 
E{n 4 Ik } — E{nQ k } = 3cr^ [3, p. 147]; and use the indepen- 
dence of nj k and ng fc , and one will obtain 

^{^4 2 } — % G n + 8^(5^ + $Q fc ) + (sj k + SQ k Y 

x 8a 4 + 8m 2 {7f<7 2 -f (C-6) 

The variance of 14 becomes 

a-l k - 4 rn 2 C 2 cr 2 + 4<t 4 (C-7) 

Finally, the variance of m $s can be written as 

= dr [4m 2 C 2 (T 2 + 4a- 4 ] (C-8) 
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Appendix D 


Evaluation of the Covariance of m p and m ss 


The covariance of m p and m 33 is defined as 


cov (m p ,m 33 ) = E{(m p — m p ){m 33 — m 33 )} 


E{m p m 33 } - m p m, 3 


As in Appendices B and C, let 


u k = 3?{n k y;j 


and 


v k = in i 2 


(D-l) 


= YaI k Y/3I k + Y aQk Y 0 Q k 
— ( s «/ t T n <*i k )( s 0i k Y 

+ (*a o. + n aQ )(s P q . + Tip ) (D-2) 


Sj t + 2 sr k ni k + n) k + Sq k + 2 SQ k riQ k + riQ k 


(D-3) 


where 


and 


Note that 


_ sine (u 0 T, s / 4) 
1 sine (u> 0 T s /2) 


sine (u> 0 T sy /2) 
O2 — — 


sine (u> 0 T s /2) 


z = u 0 T„ 




C\ = -S'Y 


1 + COS 


^ U) 0 T S y 


Then 


E{m p m ss } = 


( N N > 

t k=l 1— 1 , 


The expectation of the product of UkVk can be expressed 
as 

E{U k V k } = 

CT n + 2<T n ^ { s «J t s 0/ t + s a Qk s 0 Qk } + + s Q k ) a n 


1 

AT2 


N 


N 


N 


+ $>{£/*} E 


k - 1 


jfe = l 


AT2 




+ ^{( 5 / fc + S Q* )( S «/ fc s j3i k + s a Qfc )} 


2 M 2/^2 

» ( v 4 m C 1 


= ^ + 2cr, 


+ ^(m 2 ^) + m 2 Cf 


cos ■ 


(r-’c? 


(D-4) 


Finally, the covariance of m p and m 5J is 


cov (m p ,m 33 ) = -[E{U k V k }-m. i 


N 


[ m 2 r 2J , 41 

[m C 2 <r n + <r n 


(D-5) 

(D-6) 
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Appendix E 

Evaluation of the Partial Derivatives 


Denote the following parameters in a convenient form, 

_ 4sinc 2 (u> 0 T s /2) 

1 sine 2 (u> 0 T sy / 4) cos (u) 0 T sy /2) 

^2 ~ / rp /9 X + 2 

COS {ijJ 0 l $ y/ 1 ) 


and 


sine (w 0 T sy /A) 

1 sine (u> 0 T s / 2) 

Z — U) 0 T S y 

In terms of K\ and K 2 j the estimate SNR* defined in 
Eq. (12) is given by 


d 2 g 


dm 2 


2 A 1 A' | m p i 2K1K2 

T / 


__ (m„ - A' 2 ;7i p ) 3 (m ss - K 2 m p ) 2 

m p > tti 5 5 


2 cos 2 (z/4) 
< 7 ^ cos 2 (z/2) 


( 1 + cos l) 


2SNR ( 1 + cos -J + ^2 


(E-3) 




dm. 


Kimp 


(m sa - I< 2 r 


2 o 2 


SNR 


(E-4) 


d 2 g 


dm 2 . 


m p , m s 5 


2 K 1 

(m JS - I< 2 m p ) 3 


SNR* = 

h\rn r 
- K 2 m p 


(E-l) 


Then = SNR, and the partial derivatives of 

this estimator function are 


(3m p 


TU p > 771 $$ 


K i K 2 m v K\ 

(m ss — I\ 2 rn P ) 2 m„ - K 2 rn p 


1 

a 2 cos ( 2 / 2 ) 


SNR (l + cos 0 



(E-2) 


2^ 


SNR 


(E-5) 


d 2 g 

dm p dm ss 


2K \ K^Tip 


I< 1 


(?n SJ - /\ 2 m,0 3 (m„ - K 2 m p y 


cos(z/2)^ L 


SNR j^l H- cos 0 + 


2 ) ' C\ j 


(E-G) 


where the second lines of Eqs. (E-2) through (E-6) arc 
obtained by substituting the expressions for m p = E{jn p } 
and fn ss = E{m 3J ) from Appendices B and C. 
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